All-optical pump-and-probe detection of two-time correlations in a Fermi gas 
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We propose an all-optical scheme to probe the dynamical correlations of a strongly-interacting gas of ultracold 
atoms in an optical lattice potential. The proposed technique is based on a pump-and-probe scheme: a coherent 
light pulse is initially converted into an atomic coherence and later retrieved after a variable storage time. The 
efficiency of the proposed method to measure the two-time one-particle Green function of the gas is validated by 
numerical and analytical calculations of the expected signal for the two cases of a normal Fermi gas and a BCS 
superfluid state. Protocols to extract the superfluid gap and the full quasi-particle dispersions are discussed. 
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I. INTRODUCTION 

Many-body quantum systems exhibit truly remarkable fea- 
tures such as high-temperature superconductivity and the frac- 
tional quantum Hall effect. Traditionally, these phenomena 
are studied in the solid state. However, in recent years di- 
lute, yet strongly interacting, atomic gases have started pro- 
viding a novel class of systems to investigate this fascinating 
physics. Their outstanding cleanliness, control, and precise 
microscopic understanding will push forward the fundamen- 
tal understanding of quantum many -body physics iQl]. 

Strongly interacting atomic quantum gases are generally 
prepared by trapping atoms in vacuum in a magnetic or optical 
potential. This offers two remarkable opportunities: Firstly, a 
superb isolation from the environment opens the door to fasci- 
nating experiments out of equilibrium to investigate genuine 
quantum dynamics. Secondly, a variety of coherent optical 
processes are available to selectively probe the quantum sys- 
tem without being disturbed by a surrounding bulk medium. 
In particular, these optical detection techniques can provide 
repetitive and almost non-destructive in-situ measurements 
M>^- The combination of these two features makes ultracold 
quantum gases ideal systems to study the non-equilibrium and 
dynamic properties of isolated quantum many body systems. 
However, the experimental study of these properties requires 
the development of novel detection schemes that are sensitive 
to a wider variety of observables of the quantum gas, e.g. its 
multi-time correlation functions. 

The prime example of an atomic quantum system mim- 
icking the physics of the solid state are interacting fermionic 
atoms in artificial lattices structures, the so called optical lat- 
tices 101 ■ The preparation of strongly correlated states in an 
optical lattice will allow for an analog simulation of complex 
quantum many body Hamiltonians. Recently, evidence for the 
stabilization of a Mott-insulating phase has been obtained by 
looking at density related quantities of the gas The 
identification and characterization of more complex quantum 
phases requires, however, the measurement of time-resolved 
single-particle correlation functions, also called Green func- 



tions, of the form {tJjI. j.{t)'ijja-' .r' (t')) . Here ^j'!^.l{t) is the an- 
nihilation (creation) operator for the internal atomic state <t at 
position r and time t. The single-particle two-time correlation 
function reveals profound information about the macroscopic 
coherence and decoherence of the systems and keeps track of 
the subtle properties of quantum phases which are not density- 
ordered, e.g. the existence of quasi-particles in a strongly cor- 
related Fermi liquid. This same correlation function plays a 
even more crucial role in the case non-equilibrium situations: 
as the most celebrated example, the particular relaxation be- 
haviour of glasses is almost invisible in one-time correlations, 
while it can be followed in full detail by measuring the two- 
time ones . 

Up to now, the single particle equal time correlation func- 
tion out of equilibrium was investigated for bosonic atoms 
ifioll . For fermions, only the energy resolved correlation 
function of an equilibrium state has so far been probed 
by momentum-integrated lfTll - [l6ll and momentum resolved 
iflTiflsIl radio-frequency or two-photon spectroscopy. A more 
elaborated scheme for the detection of the two-time correla- 
tion function based on the immersion of a ion into quantum 
gases has been proposed in lfl9tl . 

Here, we propose an all-optical pump-and-probe scheme 
to extract quantitative information on the microscopic physics 
of a Fermi gas and in particular on its two-time correlation 
functions. A pump sequence firstly brings the system into a 
quantum superposition of its initial state and an excited state. 
The response of the system to a second probe pulse sequence 
is then measured after a variable time delay. In this way, in- 
formation on the time evolution of the atomic two-time corre- 
lations is converted into easily detectable observables, such as 
the intensity and the phase of the outgoing light. 

From an alternative point of view, our scheme can be seen 
as an application of light storage techniques ll20l - {22ll to the di- 
agnostic of many-body systems: a coherent pulse of light is 
stored in a quantum gas and retrieved at a later time after a 
variable interval. Information on the system is extracted from 
the properties of the retrieved light. Differently from standard 
light storage experiments where it is a purely detrimental ef- 
fect, decoherence of the stored pulse as a function of storage 
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time is in our scheme the crucial tool to obtain information on 
the many-body dynamics of the underlying quantum gas. 

The structure of the paper is the following. In SecHIl the 
measurement schemes are introduced and analytical expres- 
sion relating the observed signal to the many-body observ- 
ables are given. In Sec|III]an application to a fermionic system 
is discussed in detail and experimental protocols to extract the 
superfluid gap and the full quasi-particle dispersions of a BCS 
superfluid are outlined. Conclusions and future perspectives 
are given in SeclIVI 



II. THE MEASUREMENT PROCEDURE 



(a) 




The gas consists of a mixture of atoms in two hyperfine 
ground states g and g' that feel an identical confinement po- 
tential. Our all-optical probing scheme involves three atomic 
levels (g, e, and m) arranged in a A scheme as schematically 
shown in Fig.[T] The m state is a long-lived electronic ground 
state whereas the e state is an electronically excited state. With 
a suitable choice of polarization and frequency, the g' atoms 
experience a negligible coupling to the pump and probe light 
fields. 

The diagnostic scheme (Fig. |2| starts with the creation of 
a coherent excitation by adiabatically switching on a laser 
of (spatially uniform) Rabi frequency Vl™ then a weaker 
collinear laser of (spatially dependent) Rabi frequency VIi{y) 
1E5I1 . The two beams are then suddenly and simultaneously 
switched-off. The frequency wi 2 of the beams are chosen 
to be resonant with the g e and m e transitions, re- 
spectively. To ensure adiabaticity of the preparation stage, the 
switch-on of the two lasers has to be performed on a time- 
scale long as compared to the internal atomic dynamics and to 
the Rabi frequencies fii, fi™. On the other hand, the switch- 
off has to be much faster than all these frequencies. 

Provided the whole excitation stage takes place on a fast 
time scale as compared to the many-body dynamics, the 
atomic position can be considered as fixed and the optical pro- 
cess described by a single atom picture where the atomic \g) 
state is adiabatically transformed into a dark state ll26l - l29ll 



FIG. 1 : (Color online) Diagram of the internal atomic levels involved 
in the proposed detection scheme. 



(b) 




FTG. 2: (Color online) Snapshots of the measurement procedure. 
From left to right: adiabatic storage of coherence from incident 
beams, free many-body evolution, light re-emission and detection. 

Even though the application scope of our measurement 
procedure is much wider, we focus our attention here onto 
the case of fermionic atoms in a two-dimensional lattice ge- 
ometry. This geometry lies at the heart of quantum simu- 
lations with the aim of exploring the mechanisms underly- 
ing high-temperature superconductivity |23]. A tight optical 
confinement potential freezes the atomic motion into a sin- 
gle xy plane. Additionally, a periodic optical lattice poten- 
tial is applied along the x andy directions to generate a two- 
dimensional lattice structure ll24ll . 
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which is fully decoupled from the excitation lasers. All other 
bright eigenstates are energetically separated and do not get 
mixed with the dark state provided the switch-on phase is per- 
formed in a slow enough manner Assuming that the phase of 
the ill 2 Rabi frequencies (once the carriers at a;i.2 are elimi- 
nated) is constant during the whole sequence, the global Berry 
phase 



I I dt (dark(i)| — Idarkft)) 
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acquired by the atom during the adiabatic evolution is easily 
shown to vanish. 

In the I ill I ^ l^^"! limit under investigation here, the ef- 
fect of the excitation stage on the initial many -body state |(/)o) 
in the lattice representation can be expressed as the transfor- 
mation: 



/ J /-vin Tm 



.)l0o). 



(3) 



Here, the V'o-.r; (V'J-r ) lattice operators destroy (create) a 
fermionic atom in the a = g,m state at the lattice site r;, 
respectively, ili.2 are the Rabi frequencies in the lattice rep- 
resentation. Initially, only the g and g' states are assumed to 
be occupied. 

The use of a locally focused laser of amplitude fli(r) on 
the (7 — > e transition allows one to selectively address a well 
defined region of the sample. The space selection is useful to 
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eliminate inhomogeneous broadening effects due to the spa- 
tially varying density in e.g. trapped systems Isollsill . If no 
spatial selection is performed, the final signal would include 
the contributions of different regions of the system. 

Once the two beams are switched off, the system evolves 
according to its many-body Hamiltonian for a storage time ts 
from the state \(j)d) to the new many-body state 

where Vmbif-s) is the many body time-evolution operator and 
the system Hamiltonian involves three contributions: is 
the Hamiltonian acting on the states g and g', the Hamil- 
tonian of the atoms in the m state, and i/int contains the in- 
teraction processes between the g, g' and m states. This time- 
evolution will change the coherence between g and m present 
in the prepared dark state. The g, m coherence remaining at 
the end of the storage time is finally probed. 

The detection of the remaining coherence can be achieved 
by different schemes: Either a fast tt pulse is applied to coher- 
ently transfer all atoms from the m to the e state and then 
coherent photons emitted on the e — 5- y transition are de- 
tected. Or the excitation is slowly released by means of a 
weak field of frequency cj2 and Rabi frequency in the contin- 
uum rij"* <C 7e that transfers the atoms adiabatically from 
the m state into a coherent superposition of m and e. In both 
cases, the electric dipole that is responsible for the emission 
at frequency wi on the e ^ g transition is proportional to the 
coherence between the g and m atomic states, 

d(r) =i5l^t^r^™,r. (5) 

The constant D depends on the details of the process and de- 
termines the duration in time of the released pulse = 
7e |(ieg/I?p, where 7e is the radiative decay rate of e state 
atoms. In the tt pulse case, the constant D is equal to the elec- 
tric dipole matrix element between the state g and e, D = dge ■ 
In the case of a slow release, D is approximately given by 
D = 2ifl2^^*^dge/je- In order for the many -body dynamics not 
to interfere with the release process, the time duration of 
this latter has to be shorter than the characteristic time scales 
of the many-body dynamics. 

The near field pattern of the emitted light amplitude is de- 
termined by the expectation value of the local dipole operator 
on the final state \(l>{ts)) 

d{r,ts) =ms)\dirmts)) (6) 

We switch now to the lattice representation by relating the 
field operators in the continuum to the lattice operators via 
the Wannier functions Wa (r) as following t/ja-^r ~ J2i '^<y l'' ^ 
ri)'/''T,ri- We checked numerically that for tight atomic Wan- 
nier functions on deep lattices, we can keep only the terms 
with Wannier factors taken at the same sites and neglect all 
other contributions. This leads to the following expression for 
the dipole operator 

d(r,t,) ~Z?^VFr.(r)(0(i,)|V;t r^i„,rJ0(is)), (7) 



with Wri(r) = w*g{Y~Ti)Wm{v—Ti) . Inserting the expression 
(HJi of the final state and switching to the Heisenberg repre- 
sentation for the operators V^^j^r^ = Ul^^{ts)lpx,riUmb{ts), 
this has the form 

d(r,t,)= Dj2Wr,{r) 

i 

= D ^Tyr,(r)(0<i|V;j,r.(i.)i»,r.(ts)|</)d>. (8) 

i 

Inserting into (O the explicit expression (O for the dark 
state and taking into account that no atoms were initially 
present in the m state, this expression can be written in the 
compact form 

that only involves a time-dependent correlation function taken 
on the initial many -body state |0o)- 

As no atoms are initially present in the m state, the initial 
many -body state \(po) exactly factorizes in a complex many- 
body state for the g, g' subspace and vacuum for the m one. 
Assuming that the few atoms that are transferred into the m 
state during the preparation stage do not significantly interact 
with the majority of atoms left in the g, g' states ifTsi [32l[33ll 
allows us to neglect the m — g,g' interaction term of the 
Hamiltonian Hint in (4) and write the time-evolution operator 
U,nb in the factorized form U,nbits) = e-^Hot^/h^-iH^tjh^ 
As a direct consequence, the Heisenberg evolution of the 
V'm.r; (t) operator is only determined by the Hm part of the 
evolution operator, while the ipg.nit) and ipg'^nit) operators 
evolve with the many-body Hq Hamiltonian in the g, g' space. 
These simple facts allow one to rewrite the dipole expectation 
value in the final form: 

x(0o|V^Ir.(i.)4,r,(O)|0o). (10) 

where all the expectation values are to be evaluated on the 
initial many-body state before the preparation stage, with no 
occupation in the m state. In particular, the m state propagator 
describes the free-particle evolution in the lattice potential. 

The far-field pattern in a direction 9 is proportional to the 
spatial Fourier transform of c?(r, tg) evaluated at a wavevector 
k equal to the projection of the emission wavevector 9 uji /c 
along the xy plane. Here, c is the velocity of light. This leads 
to the following expression for the far-field emission ampli- 
tude at a distance 11 = R 9: 

q 

(11) 
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where Wm(q) is the free-particle dispersion of m state atoms 
in the lattice potential and the coefficient Ck is defined as 

_ Dujfn,{k)w{k) 

Invariance under translations along the plane guarantees that 
the coherent emission amplitude in the 6 direction, i.e. 
with an in-plane wavevector k, only depends on the inci- 
dent probe amplitude rii(k) at the same k. Here we have 
set r2i(k) = J2j f^i(r)e^*'^ '"3, we have defined i/'g q = 
(7V)-i/2 ^^^^^g-«q and we have used a lattice with 
N sites neglecting boundary effects. The factor VF(k) = 
J (Prw*{r)'Wm{'t^)e~^^''^ is a slowly varying envelop stem- 
ming from the tight atomic Wannier functions. 

Expression ( fTTI ) relates the coherent amplitude of 
the released light to the time-dependent one-body Green func- 
tion of a generic many-body gas. It is one key result of the 
present paper. In the limiting case oj„i (q) = uj°^ where the m 
atoms do not appreciably move during the time tg, the far- field 
amplitude ( fTTI ) can be further simplified into the form 

= Cke-*""*=(0o|V^Ir.(^.)4,r.(O)|0o), (13) 

which only involves the local value of the Green function of g 
atoms. 

Experimentally, the coherent E^^^ amplitude can be mea- 
sured by homodyne detection of the emission with a stronger 
reference beam at oji. The intensity and phase of is in- 
ferred from the amplitude and phase of the oscillations in the 
interference signal as a function of the mixing phase. This 
procedure requires coherence at the g — > rn frequency which 
can be easily achieved if all f^i, ilj" °"* fields are obtained 
from a single laser source. 

Another quantity of interest is the total (i.e. coherent and 
incoherent) intensity pattern in either the far- or the near-field. 
Differently from the coherent amplitude (fTTT i. these involve 
higher order correlations of the many-body gas. For instance, 
the near-field dipole pattern /(r) = ((i^(r) d{r)) reads: 

X (vac| Vi„,r, (0) 1^;^^r. {ts)4'm,r, (ts) ll>L.r, (0)|vac) 
X (0o|^lr, (0) 4,r. its) ^Ir. V^ff.r, (O)|0o). 

For a localized beam fli [rj], the /(r) signal is proportional to 
a fixed envelope determined by the motion of atoms in the m 
state times a two-body Green function of g atoms. The corre- 
lation function of the g state can be understood as measuring 
the density at time tg at site Ti if one has removed an atom at 
time at site rj. This scheme looks promising e.g. to follow 
the dynamics of holes in anti-ferromagnetic states ll35i[36ll . 



for a weakly attractive, unpolarized two-component Fermi gas 
in an optical lattice at half filling. In particular we show 
how the proposed method is able to identify a superfluid state 
and its quasiparticles from the measured two-time correlation 
function. 

In the normal state, the dispersion relation of quasiparticles 
is given by the free-particle dispersion in the lattice. Here 
we take the tight-binding form ftwg,g'.m(q) = g' m ^ 
2 Jg,g',m[cos((7^a) + cos((7j,a)]. While the g, g' atoms feel the 
same potential, Jg ~ Jgi, the hopping J,„ for the m state 
atoms can be different. In the following we set oj", ^ uj° — Q 
and focus on the case of half-filling. 

In the superfluid state, the quasiparticle dispersion pre- 
dicted by BCS theory consists of two branches E'^ = 

± •y/[^a^^"(q)p"+"A^ separated by a gap of amplitude 2 A. The 
one-body Green function 0g(q, t) ~ (ip^^ ^{t) tl^g.qiO)) for the 
BCS phase reads Q 

g,(q,i)=^4/(£;+)e'("-+^^/'')* 

The Bogoliubov coefficients are defined as Uq, Wq = 
i [l + h(jjg{q)/ E^'\ and the Fermi distribution as f{E) = 
(1 + e^/'=«'^)"^ is the mean-field shift iH. In what 
follows, we shall focus our attention on low temperature T 
for which the upper branch E^ is almost empty and can be 
neglected. Under such an assumption the emission amplitude 
(fTTI) becomes 

q 

Its Fourier transform with respect to the storage time tg has 
the form 

^r^) = ^^t'^/(£;q)<5K-^„,(q+k)+C..„f + i?q/n). 

q 

For each value ujs of the frequency, the signal comes from the 
wavevectors q which fulfill 

LUg ^ - w,„f + r ujg (q + k) - E' /h. (14) 

In the following we will neglect the contributions by "Msf 
since these can be eliminated in the homodyne detection HUR . 
Several regimes can be identified depending on the value of 
the hopping ratio r = Jm/ Jg- Experimentally, the hopping 
amplitudes can be varied within some range by tuning the fre- 
quency and polarization of the lattice beams, or, if necessary, 
by using more complex multi-photon transitions instead than 
the simple Raman scheme discussed so far. 



III. APPLICATION TO BCS SUPERFLUID 



A. Small hopping ratio r <g; 1 



In order to demonstrate the efficiency of the proposed de- 
tection technique, we now calculate the signal that is expected 



The physics is the simplest in the r <C 1 case where the 
atoms in the m state do not move during the experiment and 
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B. Equal hopping ratio r = 1 
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FIG. 3: (Color online) Time-dependence of the emission amplitude 
in the normal (k — 0) direction for a normal state N (blue solid 
line) and a BCS superfluid with gap A /4 Jg =0.2 (red dashed line). 
Hopping ratio r = 0.1. Temperature ksT = Jg/50. 
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FIG. 4: (Color online) Frequency-dependence of the emission am- 
plitude in the normal (k — 0) direction for a normal state N (blue 
solid line) and a BCS superfluid with gap A/4 Jg = 0.2 (red dashed 
line). Hopping ratio r — 0.1. Temperature ksT = Jg/50. 



the emission amplitude is determined by the local Green func- 
tion ( fT3b . As one can see in Fig. [3] the emission amplitude 
for a superfluid state as a function of storage time tg shows 
a slowly decaying oscillation at a low frequency determined 
by the BCS gap A. On top of this slow oscillation, faster 
and quickly decaying oscillations are visible at frequencies on 
the order of the Fermi energy (i.e. the band width Jg). The 
long lasting, slow oscillations are a signature of the superfluid 
state. They disappear in a normal state where one is left with 
fast and quickly decaying oscillations. 

This physics is easily understood looking at the correspond- 
ing frequency spectra plotted in Fig. |4] In the limiting case 
r ^ 0, the spectrum recovers the density of states for quasi- 
particles. In the normal state, the spectrum has a broad shape 
extending up to hojmax — 4Jg (1 — r) and showing a singu- 
larity at = as a consequence of the perfect nesting of the 
square Fermi surface at half-filling. In the superfluid state, the 
dominant feature is the peak at hajg ~ A that limits the spec- 
trum from below and from which the BCS gap is immediately 
extracted. In this state the upper limit of the signal is shifted 

to fiWmax = -4r J„ - 
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FIG. 5: (Color online) Frequency dependence of the emission ampli- 
tude from a BCS superfluid with A/4 Jg — 0.2 at different angles, 



k — kx — ky 



0, 7r/4a, TV /2a, with hopping ratio r — 1. Inset: 



16 J2 + A2. 



magnified view of the k = curve. The arrows indicate the corre- 
sponding spectral minimum ijjmin. Temperature fcsT = Jg/50. 

The coherent emission spectra in the case of equal hopping 
amplitude r = 1 show a rich structure that strongly depends 
on the wavevector k (Fig. |5]l. Even if the physics is some- 
how more involved than in the r ^ 1 case considered in the 
previous subsection, still the observed signal can be used to 
obtain useful information on the many-body system, e.g. its 
superfluid gap. 

Let us first focus on the coherent emission in the k = 
direction. At the lower boundary a large signal is found at 
/iWniin ~ (^^/4>/c()/2 for A 4Jg (see the inset) which 
originates from quasiparticles at q = 0. The long tail that 
appears at high frequencies past hujg = A is a direct conse- 
quence of the smearing out of the Fermi surface on an energy 
scale A in the BCS state. 

The emission spectrum in the direction along the diagonal 
of the Brillouin zone k = (/c, k) with k ~ 7r/2a (a is the 
lattice constant) is characterized by two peaks and a broad 
background with quite sharp edges: most visible is the strong 
peak at hojs = A that originates from the divergence of the 
density of states at the Fermi level in a BCS state. This peak 
persists for different values of k (cf. Fig.|5]A: = n/Aa) and its 
position can be used to experimentally measure the amplitude 
A of the gap. 



C. Higli Iiopping ratio r 1 

We conclude our study with a brief account of the case of 
a high hopping ratio r ^ 1. Examples of spectra for r = 3 
are plotted in Fig. |7]i: in particular, they show a clear peak 
at hujg = A independent of the direction of the light. This 
distinctive feature allows for a direct measurement of the gap 
amplitude A. 

Furthermore, the full dispersion of the BCS quasiparticles 
can be extracted from the position of the lower edge of 
the spectrum. For r 3> 1 the r-dependent term in Eq. ( fT4l ) 
dominates and determines the q values that correspond to 
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IV. CONCLUSIONS AND PERSPECTIVES 
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FIG. 6: (Color online) k-dependence of the lower edge of the spec- 
trum compared to the quasiparticle dispersion of the BCS su- 
perfluid with A/4Jg = 0.2. From top to bottom, hopping ratio 
r = 1, 1.5, 3. Curves for different r are offset by 0.5 for better visi- 
bility. 
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FIG. 7: (Color online) Frequency dependence of the emission ampli- 
tude from a BCS superfluid with A/4 Jg — 0.2 at different angles, 
k = kx = ky = Q, 7r/4a, -k /2a. Hopping ratio r = 3. The ar- 
rows indicate the corresponding spectral minimum ixJmin. Tempera- 
ture kBT = Jg /50. 



In summary, we have proposed a novel all-optical, spatially 
selective and almost non-destructive technique to probe in situ 
the microscopic many-body dynamics of a gas of interacting 
ultracold atoms. The technique is inspired to recent light stor- 
age experiments and is based on the creation of an atomic co- 
herence by coherent absorption of a pump laser pulse and its 
later retrieval after a variable storage time: information on the 
many -body dynamics of the quantum gas is extracted from the 
amplitude and coherence properties of the retrieved light. Dif- 
ferently from most previous measurement schemes, the use of 
a spatially localized pump spot will allow to individually ad- 
dress the different coexisting quantum phases that can appear 
in a trapped system. 

The efficiency of the proposed measurement scheme is 
tested on the specific, analytically tractable example of a two- 
dimensional BCS superfluid. Protocols to extract the su- 
perfluid gap and the quasi-particle dispersion are presented, 
which take into account some most significant difficulties that 
arise from the internal structure of the atoms. 

As our scheme consists of the measurement of two-times 
correlation functions, it is expected to be of great utility in the 
study of the non-equilibrium dynamics of a quantum system: 
on one hand, its almost non-destructive nature suggests that a 
series of many measurements can be performed at a high repe- 
tition rate without significantly perturbing the system dynam- 
ics. On the other hand, the observed quantities play a crucial 
role in the characterization of relaxation dynamics [ 9, 381: for 
instance,they may serve to identify the glassiness of a system 
in the presence of disorder 

Future work will investigate the extension of the method 
to more complex, three-dimensional geometries: differently 
from the two-dimensional geometry considered so far, this re- 
quires a careful treatment of light propagation across a bulk 
sample in both the excitation and the retrieval stages. Prelim- 
inary work in this direction has appeared as IjtIi . 



the spectral edges: the contribution of quasi-particles with 
momentum q = — k determines the sharp lower edge at 
^min = -4rJg - EZ]^- 

The dependence of the lower spectral edge on the emission 
direction {k = = ky) is shown in Fig.|6]for different values 
of r and compared to the quasi-particle dispersion. While the 
agreement is limited to the special points k = 0, 7r/2a for 
r = 1, it quickly improves for larger r; a reasonably accurate 
image of the quasi-particle dispersion around fc = 0, 7r/2a is 
already recovered for r > 3. 
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